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Abstract 

Phenomenology of neutralino dark matter in the minimal supersymmetric model is discussed 
for a scenario where the lightest Higgs boson mass is lighter than 114.4 GeV. We show that the 
scenario is consistent not only with many collider experiments but also with the observed relic 
abundance of dark matter. The allowed region may be probed by experiments of B s — * u + [i~ in 
near future. The scenario predicts a large scattering cross section between the dark matter and 
ordinary matter and thus it may be tested in present direct detection experiments of dark matter. 

PACS numbers: 95.35.+d, 12.60.Jv, 14.80.Cp, 13.66.Fg 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) experiments will start soon. The primary goal of the 
LHC experiments is the search for the Higgs boson 0, which is the last undiscovered particle 
in the standard model of the particle physics (SM). The Higgs boson is the only scalar 
particle in the SM and rules masses of the other particles via couplings with them. Hence, 
the discovery of the Higgs particle is indispensable to understand the nature. 

The present bound on the Higgs mass in the SM comes from the LEP2 experiments and 
is given as [l| 

m lj> > 114.4 GeV (95%C.L.). (1) 

In addition to the bound, the Higgs mass should be within TeV scale theoretically because of 
the partial wave unitarity of weak gauge boson scatterings. On the other hand, once quantum 
correction to the Higgs mass square is taken into account, it has quadratic divergence. 
Hence, the natural scale of the SM Higgs mass is the cutoff scale of the SM, which is, for 
example, the grand unification theory (GUT) scale {Ma ~ 10 16 GeV) or the Planck scale 
(~ 10 19 GeV). This problem in the SM is called as the hierarchy problem and considered to 
be very serious. To solve this problem has been the paradigm for particle phenomenology 
for a couple of decades. 

The supersymmetry (SUSY) [2j is one of solutions to this problem. The quadratic di- 
vergence is canceled between loop diagrams of SM particles and those of SUSY partners. 
Moreover, SUSY theory also has some advantages; the lightest SUSY particle (LSP) is pro- 
vided as a candidate for dark matter with R-parity conservation, and the grand unification is 
improved in comparison with the SM. By these reasons, SUSY theory is the most celebrated 
candidate for physics beyond the SM. 

In the minimal SUSY SM (MSSM), two Higgs doublet fields (H u ,Hd) are required to 
provide mass terms for quarks and leptons. Hence, the MSSM predicts three neutral and 
two charged Higgs bosons: two CP-even scalar ones h (lighter), H (heavier), one CP-odd 
pseudoscalar one A, and a pair of charged ones H ± . The mass of the lightest one, m^, can 
be lighter than the SM bound ([1]) if gzzh, which is the coupling of h to the Z boson, is 
significantly smaller than the SM value si, 0, S]. We call the MSSM with the Higgs boson 
whose mass is lighter than 114.4 GeV as the light Higgs boson scenario (LHS). The LHS 
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might explain the 2.3 a level excess of events at 98 GeV in the LEP2 experiments 
though it is very difficult to be explained in the SM with = 98 GeV. 

A reference model widely used to study the MSSM is the constrained MSSM (CMSSM), 
which is parametrized by mo, miw (universal scalar and gaugino masses at Mq), A q (uni- 
versal trilinear coupling at Mq), tan/3 = (H u )/(Ha), and sign(/x). Unfortunately, the LHS 
cannot be realized in the CMSSM. Thus, studies for the LHS has been restricted in spite 
of its potential importance 4, . In recent works of the LHS, phenomenological 

aspect with m h ~ 98 GeV 6| and m h < 90 GeV [7,] and a solution to the little hierarchy 
problem |8| are discussed. 

In this Letter, we investigate phenomenology of neutralino dark matter in the LHS. In 
particular, it is studied whether the LSP in the LHS is a viable dark matter candidate, 
that is, whether the parameter region where the relic abudance of the LSP is consistent 
with the observed abundance by the WMAP experiment 10j is allowed by the recent severe 
constraints from b — > 57 and B s — > We calculate the relic abundance by using the 

SUSY model with non- universal scalar masses for the Higgs multiplets (NUHM), which is 
the simplest example realizing the LHS and consistent with GUT. We also calculate the 
cross section of the dark matter with ordinary matter and compare with present bounds by 
direct detection experiments of dark matter. 



II. LIGHT HIGGS BOSON SCENARIO 

The mass eigenstates for the CP-even neutral Higgs bosons in the MSSM are given by 



— sin a cos a 
cos a sin a 




(2) 



where a is the angle diagonalizing the mass square matrix of the Higgs bosons, 



m\s 2 p + m 2 z cj + A dd 
-(m 2 A + m 2 z )s p Cf3 + A 



du 



-{m 2 A + m 2 z )s p cp + A du 



(3) 



Here, m A {m z ) is the pseudoscalar Higgs boson (Z boson) mass, c@ = cos/3, sp = sin/3, and 
Ajj denotes the radiative correction, which is given in Ref. j^j]. 

The LEP collaborations have analyzed their data using several benchmark scenarios to 



search Higgs bosons in the MSSM 



3j]. In their results, it is implied that < 114.4 GeV is 
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allowed if gzzh is significantly suppressed since the primary mode to search h in the LEP2 
experiments is e + e" — > Zh. The coupling gzzh is given as gzzh = gzz<t>sm(/3 — a), where 
gzz<t> is the coupling of to Z. To avoid the constraint dTJ, so small sin(/3 — a) is not needed. 
When 90 GeV < < 114.4 GeV, the Higgs boson h with sin(/3 — a) < 0.5 can evade Higgs 
boson searches by the LEP2 experiments. For this suppression, a large mixing a between 
h and H is required to cancel /3, and hence m# ~ m^. Note that m# > 114.4 GeV is also 
required because H has not been discovered yet. On the other hand, suppressed sin(/3 — a) 
leads unsuppressed cos(/3— a). This results in large gzAh, since #za/i oc cos(/5— a). Therefore, 
we should be also careful in e + e~ —* Ah mode in the LHS. However, since the gzAh coupling 
originates in a derivative coupling, no additional constraint appears for m h > 90 GeV due 
to P-wave suppression {j}]. 

The LHS cannot be realized in the CMSSM as mentioned above. If mo and mi/2 are 
fixed to derive mh < 114.4 GeV, masses of other SUSY particles are too small. For re- 
alization of the LHS, the masses of the Higgs doublets should be different from others as 



m 



H u (H d ) = (1 + ^-fMfld)) m o- ft * s reasonable because the Higgs multiplets are not necessarily 
in the same multiplet of GUT with other scalar particles. The simplest model with this 
boundary condition is the NUHM. Hence, we adopt this model as a reference model to in- 
vestigate the LHS. The NUHM has six free parameters, (mo, mi/2, Aq, tan/3, //, m^), where 
(mo, mi/2, A ) are defined at M G and the other is defined at the electroweak scale. This 
parametrization allows two Higgs doublets to take their masses at Mq as free parameters. 
Using these values at rriz scale, a boundary condition at Mq is derived by the renormal- 
ization group running. Then, masses of SUSY particles run back from Mq to mz- In this 
work, the renormalization group running is evaluated by IS A JET 7.75 11]. 



III. CONSTRAINTS 

We consider the following constraints from collider experiments. 

(i) LEP2 Zh/ZH constraints: In order to satisfy this constraints, sin(/3 — a) < 0.5 is 
required for 90 GeV < nih < 114.4 GeV. In our analysis, this constraint is imposed in a 
more precise form dependent on m^ with assuming that the branching ratios of h decays 
are the same as those of the SM This assumption is good enough because SUSY loop 
corrections to the Yukawa coupling of b quark are small in our analysis in which we take 
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tan/3 = 10 in order to satisfy the constraint from B s — > fi + fi~ process (See (v) also). In 
addition, m# > 114.4 GeV is required though m# cannot be very different from rrih- Even for 
rrth < 90 GeV, the Higgs boson can evade the LEP2 searches with smaller sin(/? — a) ~ 0.2. 
However, this region is very restricted, and therefore it is ignored in this work for simplicity. 

(ii) LEP2 Ah /AH constraints: In our calculation, the constraint on gzAh depends on 
the masses itia and rrih |3j. We assumed m^-max benchmark for the final states of the 
Higgs bosons decay. Since the dominant decay mode in our case is b quark pair production, 
which is the same as the m^-max benchmark, the uncertainty from this assumption is small 
enough. For 90 GeV < rrih < 114.4 GeV, this constraint is not severe due to the P-wave 
suppression. 

(iii) SUSY particle search constraints 1121 : We consider the constraints for masses of 

n 

SUSY particles through micrOMEGAs 1.3.7 package [13J which is also used to calculate 
the following two processes and the abundance of the LSP in the LHS. In particular, 
the constraint for the lightest chargino mass is precisely given depending on the lightest 
neutralino and sneutrino masses. 

(iv) b — > S7 process: For the small charged Higgs boson mass rriH±, the constraint from 
Br (b — > 57) is serious, since H ± contribute constructively to the SM prediction [ijj]. The 
charged Higgs boson mass is predicted as the electroweak scale in the LHS. This is because 
generally m H ± ~ mu ~ tua in the MSSM and rrin ~ iTih is required for large mixing. 
Hence, the contribution from H ± should somewhat be canceled by SUSY contribution in 
the LHS 8j. In particular, the H ± contribution can be compensated by large A-terms 1 . 
The present experimental bound is reported as [lf| 

Br(6 -> 57) = (355 ± 24+? ± 3) x 10~ 6 . (4) 

In our analysis 2 , we search allowed regions for 304 x 10 -6 < Br(6 — > S7) < 406 x 10 -6 . 



1 Large A-terms may cause the problem of charge (and color) breaking minima 15[. We checked that 
values of Aq used in our analysis are acceptable because the MSSM minimum is stable enough even if it 
is not the global one. 

2 The recent theoretical prediction of Br(6 — > s-f) in the SM is significantly lowered as Br(6 — > s-f) — 
(2.98 ±0.26) x 10 -4 17j . This improvement was not included in our calculation which gives Br(b — > s-f) — 



3.62 x 10 -4 as the SM prediction. Since the deviation of the new value from the present experimental result 
(TJ| is larger than that of the value we used, it will shift the allowed region we obtained for Br (b — > sj) to 
a heavier mass region in (7711/2, m o) plane (See also Fig.|TJ) and/or enable us to use smaller tan/3. 
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(v) B s — > n + \i~ process: The process, B s — > is mediated by neutral Higgs bosons. 

The H/A contribution to B s — > is proportional to ra.^ 4 . In the LHS, this contribution 

is large since ra^ ~ ra# ~ ra^. Moreover, this contribution is proportional to (tan/?) 6 . As a 
result, the LHS with large tan/? > 20 is almost forbidden. The present experimental bound 
is given by [3] 

Br(B s -> fi+ij-) < 5.8 x 10~ 8 (95 % C.L.). (5) 
This upper bound is used for our calculation. 

IV. RESULTS FOR THE DARK MATTER 

The purpose of this work is not to investigate all parameter space exhaustively but to 



indicate that t 
experiment 



re LHS is consistent with dark matter abundance observed by the WMAP 



10|, 



n BM h 2 = (0.1277±°;°o7g) - (0.02229 ± 0.00073), (6) 



101 ]. We adopt the allowed 



where h = 0.73^q;q 4 stands for the reduced Hubble constant 
region as 0.0895 < Odm/j 2 < 0.1214. First, we take tan/3 = 10. Small tan/3 is favored by 
the B s — ► process. On the other hand, too small values of tan /3 are disfavored by the 

b — > S7 process with moderate values of A for the stable MSSM vacuum. Furthermore, we 
choose \i = 300, 600, 700 GeV, where only positive fi is considered because it is favored by 
muon g— 2. The values of Aq are chosen so as to cancel the contribution to b — > 57. Then, 
we scan the parameter region where 50 GeV < ran, my 2 < 700 GeV and ra^ = 80-140 GeV. 
This range of ra^ is a typical one in the LHS [3, Q] ■ 

In Fig.Q], we fix ra^ = 100 GeV for simplicity. The region consistent with the dark 
matter abundance observed by the WMAP experiment is shown for \i = 600 GeV and 
A = —1000 GeV in (ra!/ 2 ,ra ) plane. The magenta (gray) region is consistent with the 
observed abundance. The strip at lower-right portion is the stau coannihilation region and 
the vertical strip around rai/2 = 150 GeV is the s-channel pseudoscalar Higgs region. In 
the former region, the abundance of the LSP is decreased sufficiently by the coannihilation 
with stau because of the small difference between their masses. Note that the difference 
is about 10% of their masses and only a mild tuning is required. In the latter region, the 
LSP annihilation is governed by the pseudoscalar mediated diagram. It is different from the 
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FIG. 1: The magenta (gray) region is consistent with observed dark matter abundance for \i = 
600 GeV, A = -1000 GeV, tan/3 = 10, and m A = 100 GeV. See the text for the explanations 
about other regions and lines. 

funnel region where the LSP annihilation picks up the pseudoscalar Higgs pole, because 
is smaller than twice of the LSP mass. Due to the smallness of m^, the annihilation cross 
section of the LSP is sufficiently large even for large mo. The yellow (light gray) region is not 
favored because a charged particle is lighter than the neutralino. The brown (black) region is 
forbidden by the appearance of tachyonic particles at the electroweak scale due to large Aq. 
The prediction for Br(6 — > sj) is denoted by solid lines, (406, 304) x 10~ 6 from top to bottom. 
It is seen that the constraint from the b — > 57 process removes broad parameter space. The 
cancelation between contributions of SUSY particles and H ± constrains m and m^j from 
both above and below. This is a salient feature in contrast with other SUSY models. The 
dashed lines shows Br(B s — > /U + /i~) = (1.5,2.5,4) x 10~ 8 from top to bottom. The recent 
bound ([5]) does not constrain the magenta (gray) region. Experiments of B s — > yU + /i~ may, 
however, probe large part of the region in the future because the sensitivities are expected 
to be better than Bi(B s — > /i + /i~) = 1 x 10~ 8 . 

In Fig. [2j we show the regions allowed by Q DM h 2 , Br {b — > 57), and Br(B s — > 
/x + /i~) in (mi/2,m ) plane for rriA = 80-140 GeV. The allowed regions are obtained 
for (fj,/ GeV, A / GeV) = (300,-700) (blue (black)), (600,-1000) (magenta (gray)), 
(700, —1100) (light blue (light gray)). For small /x, the large allowed region is derived and 
the portion around m 1( / 2 ~ 300 GeV is mild bino-Higgsino mixing region. If the mixing is 
large, the LSP pair can intensely annihilate and the relic abundance is too small. On the 
other hand, the allowed region is very restricted for large /1. Figure [3] shows values of for 
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FIG. 2: The regions are allowed by the recent results for the dark matter abundance, b — > 57, and 
B s -» n+n- with (///GeV, A /GeV) = (300,-700) (blue (black)), (600,-1000) (magenta (gray)), 
(700, —1100) (light blue (light gray)). The regions are searched for tua = 80-140 GeV. 
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FIG. 3: The mass mt of the lightest CP-even scalar Higgs for the allowed regions in Fig.[2j 



the allowed regions presented in Fig.[2j It is clear that a wide range of rrih bellow 114.4 GeV 
is allowed. 

The spin-independent neutralino-nucleon elastic-scattering cross section for direct detec- 
tion experiments is large for small masses of the neutral Higgs bosons [4, 19]. In the LHS, 
because all Higgs boson masses are as light as the electroweak scale, the prediction for the 
cross section is generically large. Thus, we should compare the prediction with the bound 
from direct detection experiments. In Fig.Hl the predictions for direct detection experiments 
are shown for the same parameter sets as Fig.0 The horizontal axis is the LSP mass and 
the vertical one is the spin-independent cross section for direct detection experiments. The 
cross section is calculated for y = 0-0.2, where y parametrizes the strangeness component 
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FIG. 4: The predictions for direct detection experiments of the dark matter are shown for the same 
parameter sets as Fig.[2j The regions are obtained for y = 0-0.2. See the text for the detail. 

in nucleons, to take into account that the cross section becomes smaller for smaller y. The 
value of y is calculated as 0.2 in Ref. [^(J, while y = is also allowed within errors. The dark 
green (dark gray) and green (gray) lines indicate the bounds from the XENON 10 and the 



CDMS-II experiments 



21] , respectively. The small value of \x is disfavored even if it gives 



large allowed region in Fig.[2j For /j = 600, 700 GeV, the cross section is consistent with 
experiments. Since the allowed region with the WMAP observed dark matter abundance 
is very restricted for /j = 700 GeV in Fig.[2l moderate value of /1 is preferable in the LHS. 
The regions allowed by direct detection experiments are not far away from recent bounds 
of the experiments, and hence the LHS predicts that the dark matter may be observed by 
experiments in near future. 

V. SUMMARY 



We have investigated phenomenology of the neutralino in the LHS for 90 GeV < < 
114.4 GeV. We confirmed that the abundance of the neutralino is consistent with the 
observed dark matter abundance. The constraint from Br (b — > 57) severely restrict a 
part of the region allowed by the abundance. However, we can choose A to cancel the 



dangerous contribution from H to b — > 37. The recent constraint on Br(B s — ► 



is 



avoided for small tan /3 and the large part of the predicted value is within the sensitivity in 
future experiments. It was also shown that the allowed region is very restricted for large 
\i. In the LHS, the cross section for direct detection experiments of dark matter is large. 
In particular, small /x is disfavored by the present bound. Hence, moderate value of /x, 



600 GeV or so, is favored in the LHS. The allowed region may be judged by B s — > 
and direct detection experiments in near future. 

Acknowledgments: We thank J. Hisano and K. Tobe for valuable comments. The work 
of MS is supported in part by the Grant-in- Aid for Young Scientists. 



[1] R. Barate et al. [LEP Working Group for Higgs boson searches], Phys. Lett. B 565, 61 (2003). 

[2] See, for example, M. Drees, R. M. Godbole and P. Roy, Theory and phenomenology of spar- 

ticles, (World Scientific, 2004). 

[3] S. Schael et al. [LEP Working Group for Higgs boson searches], Eur. Phys. J. C 47, 547 (2006). 

[4] A. Bottino, F. Donate, N. Fornengo and S. Scopel, Phys. Rev. D 63, 125003 (2001). 

[5] G. L. Kane, T. T. Wang, B. D. Nelson and L. T. Wang, Phys. Rev. D 71, 035006 (2005). 

[6] M. Drees, Phys. Rev. D 71, 115006 (2005). 



[7] A. Belyaev, Q. H. Cao, D. Nomura, K. Tobe and C. P. Yuan, |arXiv:hep- ph/0609079 

[8] S. G. Kim, N. Maekawa, A. Matsuzaki, K. Sakurai, A. I. Sanda and T. Yoshikawa, Phys. Rev. 

D 74, 115016 (2006). 
[9] A. Bottino, N. Fornengo and S. Scopel, Nucl. Phys. B 608, 461 (2001). 
[10] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 170, 377 (2007). 



[11] F. E. Paige, S. D. Protopopescu, H. Baer and X. Tata, arXiv:hep-ph/0312045; 
[12] W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 

[13] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 149, 103 

(2002) ; Comput. Phys. Commun. 174, 577 (2006). 

[14] P. Gambino and M. Misiak, Nucl. Phys. B 611, 338 (2001). 

[15] A. Kusenko, P. Langacker and G. Segre, Phys. Rev. D 54, 5824 (1996), and references therein. 
[16] E. Barberio et al. [Heavy Flavor Averaging Group (HFAG) Collaboration], arXiv:0704.3575 
[hep-ex] . 

[17] T. Becher and M. Neubert, Phys. Rev. Lett. 98, 022003 (2007). 
[18] P. Mack [CDF Collaboration], |ar2Gvi071L12502J [hep-ex]. 

[19] For example, A. Bottino, F. Donato, N. Fornengo and S. Scopel, Phys. Rev. D 68, 043506 

(2003) ; A. Bottino, N. Fornengo and S. Scopel, Phys. Rev. D 67, 063519 (2003). 



10 



[20] J. Gasser, H. Leutwyler and M. E. Sainio, Phys. Lett. B 253, 252 (1991). 
[21] J. Angle et al. [XENON Collaboration]. larXiv:0706.0039 [astro-ph]: D. S. Akerib et al. [CDMS 
Collaboration], Phys. Rev. Lett. 96, 011302 (2006). 



11 



